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We investigate the dynamical properties of low dimensional systems, driven by external noise 
sources. Specifically we consider a resistively shunted Josephson junction and a one dimensional 
quantum liquid in a commensurate lattice potential, subject to 1// noise. In absence of nonlinear 
coupling, we have shown previously that these systems establish a non-equilibrium critical steady 
state [Nature Phys. 6, 806 (2010)]. Here wc use this state as the basis for a controlled rcnor- 
malization group analysis using the Keldysh path integral formulation to treat the non linearities: 
the Josephson coupling and the commensurate lattice. The analysis to first order in the coupling 
constant indicates transitions between superconducting and localized regimes that are smoothly 
connected to the respective equilibrium transitions. However at second order, the back action of 
the mode coupling on the critical state leads to renormalization of dissipation and emergence of an 
effective temperature. In the Josephson junction the temperature is paramctrically small allowing to 
observe a universal crossover between the superconducting and insulating regimes. The IV charac- 
teristics of the junction displays algebraic behavior controlled by the underlying critical state over a 
wide range. In the noisy one dimensional liquid the generated dissipation and effective temperature 
are not small as in the junction. We find a crossover between a quasi-localized regime dominated 
by dissipation and another dominated by temperature. However since in the thermal regime the 
thcrmalization rate is parametrically small, signatures of the non-equilibrium critical state can be 
seen in transient dynamics. 

PACS numbers: 05.70.Ln,37.10. Jk,71.10.Pm,03.75.Kk 

I. INTRODUCTION 

Advances in the fabrication and coherent manipulation of quantum many-body systems, such as mesoscopic 
devices^ , superconducting cireuits^"^ and ultracold atomic gases^ allow to investigate new regimes of many- 
body physics and open the way to interesting technological applications. A recurring issue with these systems 
is that they are often being driven out of equilibrium, either purposely as part of the preparation scheme, or 
due to their sensitivity to external noise soiirces. Such studies provide a strong motivation for understanding 
the properties of quantum many-body systems under non-equilibrium conditions (sec e.g. [6,7]) . Since a full 
microscopic description of such complex systems is in general hopeless, it is important to identify regimes 
where the complexity of the dynamics gives rise to robust emergent behavior that is independent on the 
microscopic details. 

In this paper we describe universal quantum phenomena that occur in many-body systems subject to 
external noise sources that are not in thermal contact with the system. We have previously identified a class 
of quantum-critical steady states that arise under these conditions**. Specifically, we considered quadratic 
model systems, such as quantum RC circuits and Luttinger liquids at zero temperature subject to 1// 
noise. Clearly the correlation and response of the systems are modified by the noise and show manifestly 
non-equilibrium phenomena. But, as we showed, the scale invariance and hence the algebraic decay of the 
correlations is not destroyed. Non linear perturbations such as a Josephson junction in the RC circuit or 
a commensurate periodic potential in the Luttinger liquid were considered using simple scaling arguments. 
This analysis showed that the perturbations can be tuned from being relevant to irrelevant by varying the 
noise power, suggesting a possible non-equilibrium phase transition. The key issue that was not addressed 
by the simple scaling arguments is the back-action of the non-linear perturbations on the quadratic elements 
of the system. 

Here we extend our previous analysis and use the critical steady states as a starting point for a perturbative 
real-time renormalization-group (RG) approach. This method was used in Rcf. 9 to describe a quench of 
a Luttinger liquid in a commensurate periodic potential. In the quench problem^'^" the mode-coupling 
induced by the lattice potential led to the generation of finite dissipation and temperature. In the present 
paper, we generalize the RG scheme and extend it to our case of interest namely the steady state regime of 
quantum systems subjected to 1// noise. This allows us to study the effects of small perturbations around 
the "non-equilibrium quantum critical points" in a controlled way. The first order expansion agrees with 
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the predictions made in Ref. 8 by means of the scaHng argument. However, the second order shows the 
appearance of new physics, which was hidden by the trivial scahng analysis. 

We highlight the new effects in the example of a "zero dimensional" noisy Josephson junction for which 
we obtain the full dynamical phase diagram. For such a system temperature and friction are generated by 
the combination of the out of equilibrium noise and the mode coupling terms. This latter effect being absent 
in the; equilibrium system. We also analyze by the present method the noisy one dimensional system of 
interacting particles. Quite interestingly in the later case, although the microscopic starting point is quite 
different from the case of quenched systems, the initial steps of the flow are similar to Ref. 9,10. 

The plan of the paper is as follows: the paper starts in Sec. II with a definition of the models we consider, 
namely the noisy Josephson junction and the noisy ID chain. The next section Sec. Ill offers a self-contained 
description of our general RG framework. This section can be safely skipped by the reader mostly interested 
by the results and not the actual method used to obtain them. Such resiilts arc discussed in the next two 
sections. Sec. IV discusses the case of the noisy Josephson junction. In particular the dynamical phase 
diagram and its physical implications have been regrouped in Sec. IV C and IV D. In a similar way, the 
noisy ID chain is examined in Sec. V, with the physical results in Sec. VB. Sec. VI closes the article with a 
summary of our results and an outline of future research directions. 



Below we describe the specific models of the noisy Josephson junction and the Luttinger liquid considered 
in this paper. Both systems consist of three key elements, (i) A quantum system of interacting particles, 
described by a many-body Hamiltonian -ffgys- (ii) A deterministic noise source, described by a stochastic time- 
dependent Hamiltonian ifnoise(i)- (iii) A dissipative quantum bath i?bath- In what follows, we will assume 
that the dissipative bath is initially prepared at equilibrium, and due to its large size, does not significantly 
deviate from this state throughout the whole experiment. Considering only noise averaged quantities, we 
then obtain a non-equilibrium steady state (NESS), in which correlation and response functions depend only 
on the time-difference. Due to the constant flow of energy through the system, from the noise source to the 
zero temperature bath, these correlations and response function explicitly violate the equilibrium fluctuation- 
dissipation theorem (FDT) and are derived in the framework of a non-equilibrium Keldysh action, to be 
described below. 



The first system we consider here is the resistively shunted Josephson junction driven by charge noise. 
This model is defined by the time-dependent Hamiltonian 



Here 6 is the phase difference across the junction and n its canonical conjugate (the charge). Ec = 2e^/C 
and Ej — hg are respectively the charging energy and the Josephson coupling of the junction. Q{t) describes 
the stochastic time-dependent fluctuations of the off-set charge of the junction (charge noise) i?bath 
models the resistor in terms of an infinite set of harmonic oscillators^^'^^. The model can also be represented 
by the electric circuit of Fig. 1, in which the charge noise is generated by a fluctuating voltage source 



Assuming a Gaussian distribution, we average over the noise and obtain an effective action* that can be 
written in the Keldysh form: 



II. MODELS 



A. Noisy shunted Josephson junction 




(1) 



Vn{t) = CQ{t). 





(2) 
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FIG. 1: Electric circuit of a resistively shunted Josephson Junction (1). (a) The junction is at equihbrium. (b) 
The junction is driven out-of-equihbrium by time-dependent charge noise, modeled by a stochastic voltage source 

Vn{t) = CQ{t) 



Here 9 and 9 are the symmetric and antisymmetric combination of the forward (9^) and backward 
paths: 6± — 9 zL 9. For an introduction to Keldysh path-integrals see for example Ref. 17. The two fields 9 
and 9 are respectively referred to as the "classical" and "quantum" fields. 
The inverse Green's function of the bath is 



Wath- I fi^coj'^t^u -2z^u;coth(^) I ^ ' 



Here R/Rq = i?((2e)^//i) and T are respectively the normalized resistance and T temperature. For 1/f 
charge noise {Ql,Qu) = (2e)^-F/|w| and 

r-i = r ^ ^ (4) 



noise ^ -47rif ^ 

For later convenience we define a scaled noise parameter 

F = 2tt{R/Rq)F (5) 

Note that the noise contributes only to the quantum-quantum part of the action: it introduces only fluctu- 
ations in the system, without any dissipation, and therefore explicitly violates the equilibrium fluctuation- 
dissipation theorem (FDT). 

We have shown in Ref. 8 that 1// noise preserves the scale invariance of the quadratic model. As a 

result physical correlations decay as power-laws in time. Of particular interest are the correlation K,{t) 
and response function 'R,{t) of the cooper pair creation operator e*^. These are easily calculated within the 
quadratic theory to be^: 



n{t-t') = ^(e^^+(*')(e*^+W-e'^-«))«e(t)sin(^7r-^^^(t-i') (6) 

An important goal of this work is to investigate the effect of the non linear coupling introduced by the 
Josephson junction on the steady state of the system. We shall study these effects in Sec. IV using a real-time 
RG within the Keldysh framework. The general RG scheme is described below in Sec. HI. 



t-t' 



2(1+F)R/Rq 



B. Noisy Luttinger liquid 



The second example we consider consists of interacting particles in one dimension, coupled to a quantum 
bath and to a source of time-dependent noise. The many-body Hamiltonian describing the system, the bath 
and the noise is: 



1 

2^ 



dx 
u 



K {uTTdx6{x))'^ + {udx(l>{x)f - g{x) cos(20) 



+ 



1 /■ dr 

H,,M-- / —f{x,t)d,<l>{x) (7) 

■K J U 
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The first two terms are the standard Luttinger hquid Hamiltonian with Luttinger parameter K and sound 
velocity u. The fields 9{x) and — (l/7r)(/)(a;) are, respectively, the phase of the bosons and their long- 
wavelength density fluctuations^^"^"'^. They satisfy the canonical commutation relation: [—^dx4'{x),9{x)\ = 
i5{x — x'). The cosine term describes an external static potential. In particular, g{x) = g u 5{x) models a 
local impurity^^, while g{x) = g models a periodic potential at commensurate filling^^. The last two terms 
of (7) describe respectively a dissipative bath and a stochastic time-dependent noise, linearly coupled to 
the long-wavelength fluctuations of the density. In Ref. 8 we have discussed specific examples, such as one 
dimensional systems of ultra-cold polar molecules and chains of trapped ions for which this model is expected 
to be relevant. 

After averaging over the noise, we obtain the following real time action: 



^ = / ^ ^ ( I) + y ax gix) cos(2^) 



(8) 



Here 



Here wc assumed a noise spectrum F{q,uj) ~ {f*{q,w)f{q,w)) = F/\uj\, implying that the noise couples 
to smooth density fluctuations but is uncorrelated at long scales. This behavior of the noise spectrum is 
consistent with recent measurements of 1/f noise on cold trapped ions^*'^^. The parameter r] describes the 
coupling strength to the bath, which is assumed to be linear, ohmic, and at zero temperature. (For the 
physical motivations of these assumptions, see the Supplementary material of Ref. 8.) For later convenience 
we also define the dimensionless noise strength 

P = (10) 

The model (8) is analogous to the zero dimensional case presented above, (2), with one important difference. 
In eq. (2) the coupling to the thermal bath R/Rq was a dimensionless parameter. Here on the other hand 
the coupling -q has units of momentum and is therefore relevant in an RG sense. To access the quantum 
critical point, wc restrict ourselves to the case of a weak coupling between the system and the bath, by 
considering the limit of ry — >■ (while keeping a fixed F). Our results are therefore valid only for length- 
scales smaller than 1/r]. In fact, this limitation is often less restrictive than other infrared cutoffs already 
existing in system, such as the finite size of the sample. 

In the absence of an external potential {g{x) = 0), the action is quadratic and one can easily compute 
any physical observable. Consider for instance the crystalline order parameter 0{x) = cos{2<p{x)), which 
measures the closeness of the system to a perfect crystal. {0{x) = 1 in a perfect crystal and 0{x) = in a 
non- interacting gas). At equilibrium, this order parameter decays algebraically signaling a quasi- long range 
ordered crystal. In the presence of the noise we have previously sliown^ that the crystalline order is still 
algebraic but with a modified exponent. Specifically, the correlation and response functions of the crystalline 
order are respectively given by 

. ^(e'*^'-''-'(e'*^'"'+e'*-"-"))^( (^_^,),!,,„_,,,. ) (II) 

n{x-x',t-t') = ^(e*'^+(^''*') (e^'^+("'*)-e*^-("'*)))«sin(7ri^)e(M(t-t')-k-^'l)^(i-i') (12) 

Here a the inter-particle distance (i.e. the natural UV cutoff). If F 7^ 0, we obtain a scale invariant state 
which explicitly violates the equilibrium fluctuation-dissipation theorem: a non-equilibrium quantum critical 
state. 

The static potential introduces a non linear mode coupling term and our goal is to understand its effect 
on the critical steady state described above. Let us flrst discuss the case of a local impurity potential, i.e. 
g{x) = g u 6{x). This can be thought of as a modified Kane-Fisher problem^^ where the Luttinger liquid 
lead is subject to local 1/f noise. Similarly to the equilibrium case, this problem can be mapped exactly to 
the Josephson junction model described above. To explicitly see the mapping, we integrate-out the Luttinger 
liquid and obtain an effective action for the field (j)^ {x = 0) on the impurity site. The quadratic part of the 
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action is then simply the inverse Green's function of this field calculated in the quadratic theory (9) 

= 0) = = ^ lim ^ J (u d,) j'^'j". ^1,^,. (13) 

Gr{uj,x = 0) = i{^l{Q)<l>^{Q)) = i ju dq , , \ . = ^ (15) 

Inverting the two by two Green's function matrix we obtain an action that is formally equivalent to the 

noisy Joscphson junction (2) with R/Rq = K and C = 0. However we note that the fields appearing in 
the impurity action have a different physical meaning. In Sec. IV E we will show how to adapt the results 
derived for the Josephson junction in order to obtain the current-voltage characteristics of the Luttinger 
liquid with a single impurity. 

The case of a commensurate periodic potential (i.e. g{x) = g constant) is different. We shall treat this 
problem using the RG approach in section V. For now we would like to point out the close similarity 
between the problem of a noisy Luttinger liquid in a commensurate potential and that of the steady state 
established after a sudden quench of the Luttinger parameter from Ki to Kf in a sine-gordon model. The 
latter problem was investigated by Mitra and one of us in Rcfs. 9,10. Without the periodic potential the 
steady state following the quench, as in our case, exhibits algebraic correlations with a modified exponent 
that depends on the strength of the quench^^. Indeed the correlations and the response in this state are fully 
described by a Keldysh action almost identical to (8) except for the substitution {Fq'^/\uj\) (K f / Ki — l)\uj\. 
When building the RG equations around the Gaussian action this difference turns out to be unimportant 
in the limit — >■ because of the relativistic dispersion of the Luttinger liquid. Since the flow will take the 
system away from this fixed point, some of the physics of these two problems could however correspond to 
two different fixed points. 



III. RG FRAMEWORK 



In both models described above, the action of the system can be written as a sum of a quadratic and a 
non-quadratic term. As noted above and in Ref. 8, the quadratic part of the action is scale invariant and 

gives rise to a critical steady state. This provides the basis for treating the non-quadratic terms using a 
perturbative renormalization group (RG) method controlled by the proximity to the scale invariant steady 
state. 

The philosophy of the RG scheme is similar to the equilibrium case. As visual we envision measuring 
the system with probes of decreasing frequency resolution A. We should be able to compute the measured 
dynamic correlations using an effective Keldysh action Sa that has A as its upper cutoff in frequency instead 
of the larger physical cutoff Aq. The effective action must be modified in such a way that all dynamical 
correlations at scales cv < A are the same as would be computed using the original action. The fiow of the 
coefficients of the different terms in the action with the decreasing frequency scale A reveals the physical 
mechanisms that play a dominant role at each scale. 

If the bare action was just the scale invariant quadratic action, then the lower energy effective action 5'a 
would of course be the same as the original action after a trivial rescaling of the coordinates and fields. 
By contrast, the non-linear term g couples modes of different frequencies. Therefore elimination of the fast 
modes will necessitate a change of the effective action to mimic the effect of the lost high frequency modes 
on the slow modes through the mode coupling term. As long as the coupling term g is small, this change 
is concomitantly small and can be computed by perturbation theory in g at each stage as we gradually 
eliminate shells of modes at decreasing frequency. 

We now wish to make an important technical note on the implementation of the cutoff within a real time 
approach. Within the imaginary time equilibrium formalism, A is often introduced as a sharp cutoff. This 
has the advantage of making a clear separation between fast and slow modes at each step of the RG. As we 
will see, this approach cannot be used for time-dependent problems because it violates causality^^. Instead, 
similar to Refs. 28,29, we implement the cutoff A by convolving the free Green's function in the quadratic 
part of the action Go(^) with a smooth cutoff function T\ = r(Ar). As we will show, a convenient choice of 
the cutoff function can significantly simplify the calculations. The challenge is how to perform the integration 
over fast modes when there is no sharp separation between fast and slow modes in frequency space. This is 
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done using a technique due to Nozieres and Gallet^^, which we adapt below to the non-equihbrium Keldysh 
framework. The derivation takes similar lines as in Refs. 9,10 by Mitra and one of us. Here we generalize the 
scheme somewhat in order to allow application to both the one-dimensional system and to the Josephson 
junction. The main results of this section are summarized by Eq. (26) and Eq. (40). These two equations 
give the general forms of the the scaling equations to first order and second order in the strength of the 
cosine potential. 



A. RG of a quadratic action 

As a warm up let us start by considering a quadratic action 

So = J d^ri d^ra ^{n)Go\n - r2)<t>{r2) (16) 

Here rj = {xi,ti) is a D = d+1 dimensional Keldysh coordinate, running on the forward and backward paths. 
For a quadratic action, the RG approach is just equivalent to rescaling the coordinates by r; — > (1 + dl)ri. 
The scaling equations of each parameter appearing in Gq^ are determined by their respective engineering 
dimension. For example, unitless parameters are scale invariant and define the critical fixed points of the 
RG flow. Energy scales, like a finite temperature, are always relevant in an RG sense and grow along the 
RG flow. 

At each step of the RG procedure, we need to effectively reduce the cutoff A by a small amount dA = dl A. 
For this purpose, we express the path- integral weight e*'^°['^i as the integral over an auxiliary field Sep: 



J V{5(j)) e^i't'-S't'){GK-dGK)-''('t>-6<t>)+i S4>{dGA)-''S4> (^lY) 



Here Ga = Go * Fa and we used a shortened notation in which the integration over the coordinates is 
substituted by products of matrices. The function dG\ implements the change of the cutoff and, for an 
infinitesimal shift dA, it equals to dAd\Gj\ = dlKd^GA- We now define a rescaled field (p' = 4> — 6<p and 
re- write the RHS of (17) as e*^!'^'''^'^!, where 

So[(l>', S(l>] = J d'^n d''r2 [<^'(Ga - dGA)- V + 54>{dG a)'^ 5(f\ (18) 

The RG procedure consists of integrating-out the "fast" degrees of freedom 5(j). If the action is quadratic, 
the resulting action for is just equivalent to the original one, up to a small shift of the cutoff, Ga — c^Ga = 
Ga-^a- Ifj on the other hand, the action contains non-quadratic terms, the integration-out of the fast degrees 
of freedom can have more important effects, which are the focus of our study. 



B. Perturbative RG of a non-quadratic action 

Consider now an action made by the sum of two terms: a quadratic term [4>\ and a small non-quadratic 
term Sg [(j)] ■ We assume that the cutoff is implemented in the quadratic part and transform the non-quadratic 

part using the definition (j) = (j)' + S(f>. Next, we expand the path integral weight e''^^"+'^''' in powers of 
Sg, integrate-out the fast degrees of freedom S<j) and re-exponentiate the resulting terms. This procedure 
generates a renormalized action for the "slow" fields (p' . Finally, we rescale the coordinates by r — >■ (1 -I- dl)r, 
to match the original cutoff, and derive closed scaling equations for the parameters of Sq and 5*^. 
In what follows we will apply this scheme to a generic homogeneous non-linear coupling of the form 

Sg[<}>]=g j d^'r F{<t>{r)) (19) 

Here F = F{(f)) is some non-quadratic function of the field </>. 



7 



1. First order in g 

The first order expansion of the Keldysh weight e*'^ = e*'^°('^)+''^f gives (1 + Sg)e^^'>^^\ Substituting 
^ = (p' + Scf) and integrating-out Scj), we obtain 

1 + i{Sg)g^ = l+igjd^r J ViS<P)F{^' + 5<^)e*W)(''«'^)"'(*^) (20) 
For simplicity, we now assume F(0) to be analytic and expand it in power-series of (/>: 



= l + iSg[cj>']+ig j d^r^"^^ c„(</)')'^-2y V{5<l>){6<l>fe 



(21) 
(22) 



= 1 + iSg[^'] + |dGA(0) j d^'r dlF{ct>') (23) 

Here wo expanded the binomial (d)' + 6(1))" , and we kept only the two lowest order terms in cJCa ^ dA. 

Finally, re-exponentiate the resulting term and rescale the coordinates by r — (1 + dl)r. We obtain that 
the renormalized action of (j)' differs from the original one by 

rfS(^) =Ddl + ^dl ^-^{<l>')Aj d'^r d^Fi^') (24) 

The first contribution of (24) corresponds to the normal scaling of the coupling, dictated by its engineering 
dimensions. The second "anomalous" contribution can lead to the generation of new couplings in the action, 
as well as to an anomalous scaling of the existing one. 

As expected, if F{(j)) is a quadratic function of 4>, d^^F is independent of the fields and the anomalous 
contributions in (24) vanishes (due to cancellation of forward and backward paths). If, on the other hand, 
F{(f)) is a periodic function of (f), the anomalous contribution can have the same functional dependence as 
the normal one. In the case of a cosine term, for instance, F{(l)) = cos((^) 9?F = — cos(0) and 



SgW] (25) 



The resulting scaling equation is 

dn r K f) . ^ 

(26) 



dl 



Where D = d+1 (I? = 1 in the Josephson junction problem and D = 2 for the Luttinger liquid). In general, 
Ga(0) = (0^)a is an increasing function of A. As a consequence, the second term in the RHS of (26) tends 
to counterbalance the first one, eventually making it irrelevant. This is the origin of both quantum phase 
transitions to be described in the next two sections. 



2. Second order in g 



The first order scaling equations do not include a back action of the interaction term on the quadratic 
terms. To account for such back action we go on to compute the flow equations to second order in the 
perturbation g: 



,iSo,A[0]+»Sg[0] 



1 



(27) 



As before, we express the field </> as cp' + 6(p and average over the "fast" degrees of freedom, 5(j). Next, we 
re-exponentiate together the first and second order, by using 



l + i(5pW-^(5|),^«exp 



(28) 
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Here the average denotes integration over the fast degrees of freedom 5(j), (...) = {...)^^ = f d{6(f))...e*^°'^^'^' '^'''l 
Hence the second order correction to the effective action in terms of 6' is 



(29) 



As usual, fohowing the integration out of the high-frequency modes we restore the cutoff to its original value 
by rescaling the coordinates r — ^ {1 + dl)r, and then rename the fields, by (j)' — ^ cj). 
Consider for instance the specific case of a homogenous cosine coupling 



g / d^r [cos(0+(r,t))-cos(<^_(r,t))] 



reie2e 



iei^(cc,t)+ie2<^(x,t) 



(30) 
(31) 



ei,e2=± 



Squaring this action, we obtain 



= E ^1^26364 j d'^n j d''r2 e^*-- (32) 



Here again r stands for the D = d + 1 space time corrdinate and we have defined $1234 = ei(/>i + €2^1 + 
e34'2 + ii4>2, where = ^(r^) = (j){xi,ti). 

The integration over the fast degrees of freedom and rescaling of the coordinates then gives 

(^1).^ E ^1^2^364(1 + dir J d'^r, J d^r2e^*-- (^1 _ ^1 A^^pl) (33) 



Using (28) we then find: 

dS^ ^ 
dl 



E ^i^2e3e4 / dn j dr2e^*--|^ [(*?234>o " 2(0')o] 



(34) 



This correction to the action includes a renormalization to the quadratic terms of the action that we can 
find by expanding the exponential factor in small fluctuations around its reference expectation value in the 
quadratic theory^^'^*. Let us define 5<^'^ = $2 _ ($2^^ and expand 



cos$ 



n=0 

00 



n=0 ^ ' 



= e-5<*'>0 (1 + $2 _ ($2^^-) ^ (3(^^4) 



(35) 



The second term in the brackets renormalizes the quadratic action. The other two terms are constant 
contributions that can be dropped (In fact these terms even cancel-out due to the sum over the eis). We can 
now substitute the definition of $1234 and plug back into Eq. (34). In doing so we note that the contribution 
of terms with ei = 63 vanishes. This is because the exponential e~2<* )o includes in this case the factor 
((^1 + 4>2Y)q 00. This leads to a restricted sum only over ei, 62 and 63, which gives 



dS^ 
dl 



A / dri dr2 E ^264 

ei,e2,e4 



ei(<Ai - (t>2) + £2^1 + £4^2 



xe, 



'0 dk 
We now explicitly perform the summation over Cj indexes and obtain 

dS(2) 



dl 



1 poo f 



cos(Gi?,A(r)) „' -|-?sm(G'fl,A(r)) 



(36) 

9GK,A(r) 



sm(Gj?,A(r)) ' - icos(Gi{,A(r)) 



OA 



dA 



dA 
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where we defined R = {ri + r2)/2 and r = ri — r2- Gk.a = ^nd Gb.a = *(</'(^)</'(0)) are the 

Keldysh and retarded components of the Green's function's in the quadratic action with a cutoff A. 

Eq. (37) amounts to a renormahzation of the quadratic action (2). The first row gives the RG transfor- 
mation of the Keldish component while the second row determines the renormahzation of the retarded and 
advanced components. These second order contributions can be rewritten as 

4 (m + ^) cos(GH,^(r))e^(«-^«-«-^(0)) (38) 

'-g^A J dr' {5{r - r') - 5{r)) + ^) sin(G^,A(r'))e^(«-^('-)-«-'^(°)) (39) 

Here we used (0) = i{4>^)- After a Fourier transform they assume the more compact form 

7^(u;) - 7^(0) \ , . 

Here to is the Fourier conjugate to r and /C(r) and 7?.(r) are respectively the correlation and response function 
of the operator e*'^ in the quadratic action So, already defined in (6) and (6) 



dl 
dl 



d( G,'{u)\_ A d d{r)\ 
dl\G^\ij) G:^\u) J 2^ \dA^ dA J\ 



IV. THE NOISY JOSEPHSON JUNCTION 



We now analyze the model of the noisy Josephson junction defined by Eq. (2). We first derive the RG 

equations of the Josephson junction valid up to second order in the Josephson coupling g. We shall then 
analyze the flow in order to obtain a dynamical phase diagram and predict the expected current- volt age 
characteristics of the junction in the different regimes. We will show that there are wide regimes in which 
the junction exhibits manifestly non-equilibrium but nevertheless universal behavior that is controlled by 
the non-equilibrium critical state of the noisy RG circuit. 



A. Derivation of the RG equations 



Here we give a detailed derivation of the RG equations valid to second order in the Josephson coupling 
constant g. Readers that are uninterested in the technical details can skip directly to the summary of the 
RG flow given by Eqs. (53) to (56). 

Gonsider first the quadratic action of the RG circuit without the Josephson coupling {g = 0). In this 
case, the scaling dimension of each term of the action is just given by their engineering dimensions. The 
dimensionlcss resistance R/Rq and the dimensionlcss noise strength F arc scale invariant and deflnc a non- 
equilibrium critical fixed point. However, if the resistor is kept at a non-vanishing temperature, then the 
temperature scales like energy and is always relevant in the RG sense, driving the system away from the 
above non-equilibrium critical state. The capacitance G is irrelevant in an RG sense and determines the 
ultraviolet cutoff of the theory: A = 1/RG. In what follows, we will neglect this term, and substitute it by 
a generic cutoff function Fa. 



1. RG flow to first order in g 

We now treat the effects of a weak Josephson coupling (gi <C A), through the perturbative RG approach 
presented in the previous section. The first order scaling equation of the coupling g is given by Eq. (26) 
derived in section III B for a generic model. 

We focus at first on the zero temperature case in which the pure RG circuit would be critical. In this 
case, the anomalous contribution to the scaling dimension can be easily computed from: 



-oo 





















'a) 







Rq AJo dw \AJ a Rq^ ' ^ ' 
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Here F is defined in (5). In the last identity we used the defining properties of the cutoff function, namely 
that r(a;) ^ 1 for a; ^ and r(a;) — )■ as a; — > oo. As expected the final result is independent of the precise 
form of the cutoff function. Substituting (41) into (26), we obtain the scaling equation of the Josephson 
coupling, valid to first order in g: 



dg 
dl 



1 



Ro 



(l + F) 



(42) 



This result agrees with the simple scaling arguments given in Ref. 8. 

Let us now reintroduce the finite temperature of the bath. As we noted already, a finite temperature is 
a relevant perturbation, leading already at tree level to the scaling equation dT/dl = T. Besides this, the 
most important effect of a finite temperature is on the first order renormalization of the coupling g. If we 
take T > 0, the integral in equation (41) is substituted by: 



2TJ 



d 1 { 

^^ = 7^ - (rA(w) + rA(-w)) cotgh 

(jJ^ Jo ^ ^ 
If, for instance, we choose a cutoff of the type = exp(— |a;|/A), we obtain: 

K,=2j,fd.e--f^cot,l.(^)=2 



T 

'a 



- 2 



(43) 



(44) 



Here ip(x) is the Digamma function. For small temperatures T <C A, Ki w 2, and we retrieve the zero 
temperature limit (41). For high temperatures, on the other hand, the integral is proportional to Ki « ^. 
The precise form of the scaling function is clearly cutoff-dependent, while its asymptotic behavior is general 
and gives Ki = 2 [l + ^). 



2. RG flow to second order in g 



As shown in Sec. Ill B, the second order contributions renormalize the quadratic part of the action, ac- 
cording to (40). To extract explicit scaling equations of the different couplings, we proceed with taking the 
low frequency expansion of (40). The zero frequency contribution of the Keldysh part is 



dG^^iw 
dT 



0) 



/ 



dt 



dG-^\t) 



dl 



-i^Io{R/Rq,F) 



Here the dimensionless factor Iq equals to 

lo = j dt 
= j dt 



d_ 
OA 
t 



A 



dt + 



OA 
dGK,A{0) 
OA 



m 
m 



= {l-Dg)Aj dtlC{t) 



(45) 

(46) 
(47) 

(48) 



where Dg = dlng/dl = 1 — R/Rq{1 + F) is the scaling dimension of the junction 5*^ in the scale invariant 
steady state of the quadratic action with 1// noise and the resistor at zero temperature. In writing (47) we 
used the fact that under these conditions the only scale is provided by the cutoff A and therefore JC{t) is a 
scaling function of At. We then integrated by parts and used (26) to obtain the result (48) . In appendix A 1 
we compute Iq for a specific choice of the cutoff (A15). Here we just note that 7o vanishes when the noise 
parameter F tends to zero and is positive for any F ^ 0. 

The zero frequency component of the Keldysh action Gj}{lo = 0) can be associated with an effective 
temperature that is generated in the RG flow. This connection is established by comparing to the low 
frequency expansion of the Keldysh component of a quadratic action at equilibrium i]uj coth(w/2r) = r]T + 
Oiyj"^), where r] is the dissipation (In our case 77 = Rq/R). Hence we have the mapping G-\w = 0) O r]T. 
The next term in the low frequency expansion of (38) is proportional to w^. Comparing the coefficient of 
this term with the low frequency expansion of the equilibrium form will not in general give the same effective 
temperature as inferred from comparison of the leading term. This reflects the fact that the system is not 
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truly at equilibrium. However in the RG sense the uP' and higher order terms are irrelevant. Therefore the 
asymptotic low frequency behavior is that of an effective equilibrium system. 

Note that the corresponding correction was computed using the bare action. As the flow proceeds one 
should in principle compute this correction with the full renormalized action, and thus in particular with a 
finite temperature. This is however very complicated, in particular given the frequency dependence of the 
generated Keldysh term. We will thus stick to the correction computed at T = 0, keeping in mind that 
at some scale the flow will be cut by the generated finite temperature. Note that if the system is really 
in equilibrium (F = 0) then /q = as it must be. This ensures that a zero temperature system will not 
generate a finite temperature in the course of renormalization. 

We turn to consider the renormalization of the off-diagonal components of the Keldysh action given by 
Eq. (39). The lowest order term in the low-frequency expansion is proportional to iw: 



dl 

The coefficient I\ is given by 



^ -= -iu) 



J dt t ^ -iw^hiP, R/Rq) (49) 



d_ (dC^\ 
OA \ dA J 



n{t) (50) 



h = A^ j dtt 

= -DgdlA^ j dt 1 1C{t) (51) 

This integral is calculated explicitly in appendix A 1. By comparing (49) with the original action (2) we 
find that Ii is associated with renormalization of the dissipation. Specifically we can rewrite Eq. (49) as 
drj/dl = /ig^/A^, where 77 = Rq/R. We note that at equilibrium (F = 0) Ii vanishes, which implies that 
then there is no renormalization of the dissipation. This is in agreement with Refs. 22,30,31. 

Renormalization of the dissipative part in the off diagonal component of the action implies a similar 
renormalization of the external noise term. Recall that the term in the Keldysh component was —2ir]\uj\ — 
4:TriF\ij\ and as we saw above it is not renormalized to second order in g. We can rewrite this term as 
—2i{r] + dri)\u)\ — 4:Tri{F — dr] / {2'7t))\uj\. In this way the renormalization of the resistor in the off diagonal term 
is complemented by the required renormalization in the diagonal term so as to maintain the same form as 
the original action (2). The price is having a renormalization of the noise term, which can be written as: 

The constant 7i is always positive, meaning that the second order contribution acts to reduce the non- 
equilibrium noise. Accordingly, the renormalization of the dissipation is accompanied by a corresponding 
negative renormalization of the temperature. 



3. RG equations and physical consequences 



1 



R_ 



1 



T 

2A~o 



F 



The above results can be summarized by the following scaling equations: 

dg 
dl 

dr 

~dl 
dR 

~d[ 
dF 
~dd 



^+:R^A^^°-^^Af^^ 



R' 



h 



RqK 
{a - l)F - 



'Al' 



(53) 
(54) 

(55) 

(56) 
(57) 



Here Aq « l/(i?C) is the bare cutoff of the junction. The unitless coefficients Iq and Ii are defined in (45) 
and (49). As shown in appendix A 1, it is possible to obtain analytic expressions for both coefficients, by an 
appropriate choice of the cutoff. 
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Let us now review the different features of the equations (53-56) which describe the RG flow of the noisy 
shunted Josephson junction (2) up to second order in g/A. The flow of the Josephson couphng term given 
by Eq. (53) is completely analogous to the equilibrium case. The first term in the brackets reflects the 
engineering dimension of the Josephson coupling, which has units of energy while the remaining terms imply 
an anomalous scaling dimension. At equilibrium and zero temperature {T = 0, F = 0) this equation reduces 
to the well known result of Refs. 32,33 with the anomalous contribution to the scaling dimension being 
~R/Rq. Both the temperature and the noise decrease the anomalous dimension, making the Josephson 
coupling less relevant. At this level the picture is identical to that given in Ref. 8, using a simple scaling 
argument. 

New physics enters in the renormalization of the circuit at second order in g/A. In particular Eq. 54 implies 
that a finite efi^ective temperature of order fif^/A is generated in the course of renormalization even if the 
temperature of the resistor is zero at the outset. Physically this can be understood as follows. In absence of a 
mode coupling term the energy imparted on the system by the noise is distributed exactly as dictated by the 
noise spectrum. For example the spectrum of fluctuations induced by 1// noise is proportional to |w|, very 
different from a thermal distribution. The effect of mode coupling is to redistribute the mode occupation 
and in particular it generates a constant contribution (i.e. frequency independent) to the distribution of 
fluctuations in the zero frequency limit. The spectral weight of fluctuations at low frequencies and hence 
the effective temperature, can be estimated using Fermi's golden rule. First it is proportional to 3^ /A which 
is related to the matrix elements for transitions from occupied modes at high energies to the low energies. 
Second, the generated effective temperature is proportional to /q, as defined in (45), which is related to the 
non-equilibrium population of high frequency modes. As required /q tends to zero at equilibrium F = 0. 
The value of the effective temperature is set by the ratio between the fluctuations in the zero frequency limit 
and the dissipation Rq /R. 

There is also a reciprocal process due to the mode coupling term, that acts to renormalize the dissipation as 
described by Eq. (55). It is worthwhile to write this equation directly in terms of the dissipation parameter 
7] = Rq/R as 

dl - A^^' 

Physically, when there is a non equilibrium population of high frequency modes they can absorb energy from 
the low frequency modes by changing their population and thereby serve as an effective dissipative bath. The 
added dissipation is proportional to the coupling between the modes (5/A)^ and to the coefficient /i, which 
is a measure of the deviation from equilibrium (7i = for F = 0). Note that, as shown in the Appendix A, 
Ii diverges for (1 + F)R/Rq < 1/2. In this regime, the above low frequency expansion becomes invalid. 
We expect that, in this regime, a non-analytic term ~ w'^ with /z < 1 is generated, which corresponds to 
sub-ohmic dissipation. The physical implications of this transition deserve further study, which is postponed 
to future work. 

The last equation (56) describes the scaling of the external noise. 1/f noise is marginal at the level of 
the quadratic action, but it becomes weakly irrelevant at second order in the mode coupling term. This is 
consistent with the system flowing toward an effective equilibrium state with a finite effective temperature. 



B. RG equations at strong coupling: Duality transformation 

So far we have treated the case of a weak Josephson coupling in a perturbative RG scheme to second 
order in the coupling g. The limit of a strong junction can be addressed by means of a well known duality 
transformation^^'^^, which is carried out in detail within the Keldysh framework in Appendix B. In this 
case the phase across the junction is essentially locked except in rare occasions when a quantum fluctuation 
induces a phase slip, which can be formally described as a tunneling event of a dual-particle across a weak 
dual junction. The current of dual particles is the rate of phase slips that gives the voltage drop on the 
original Josephson junction. 

We can thus regard the dual junction as embedded in a complete dual circuit that is identical in structure 
to the original circuit with the transformations. 

9 Qdual 

R R-^ 

F ^ {R/RqfF 

T ^ T (59) 
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The prefactor of the dual cosine, ^duab can be evaluated using W.K.B. approximation as^^ 5duai ~ 

Aexp (^—y/gJlCj . In particular, if g ^ A, then gduai is small, facilitating the perturbative RG approach in 

9duai- The RG equations are exactly the same as (53)-(56) above except with the replacements dictated by 
the duality transformation (59). 



C. Dynamical phcise diagram 

The RG flow equations derived above for the weak and strong coupling limits can serve to identify the 
different steady state regimes of the system. We will show in this section that the entire parameter space 
is divided into three distinct regimes, which can be best described as (i) non-equilibrium superconducting, 
(ii) non-equilibrium insulating and (iii) thermal metal. We shall identify the characteristic scales that define 
the different regimes and derive analytic expressions for the crossover lines which separate them. 

We start the analysis of the dynamical regimes from the weak coupling limit, that is when the (bare) 
Josephson coupling is much smaller than the cutoff Aq = l/(i?C). A rough separation into a superconducting 
and an insulating regime can be made using the first order flow of the coupling g determined by Eq. (53) 
with T = 



dg 



1-1^(1 + ^) 
-Kg 



(60) 



At this order we would identify the insulating phase as the region R/Rq < (1 + F)~^, where the Josephson 
coupling is irrelevant, and the superconducting phase as the complementary region where g flows to strong 
coupling. 

The picture becomes somewhat more involved when we consider the second order correction to the flow. 
Most importantly, Eq. (54) shows that in presence of external noise, the Josephson coupling acts to generate 
an effective temperature. There is always a scale A^ below which the temperature becomes the dominant 
term and where the junction will exhibit metallic behavior. But that scale may be so small that ^(A) 
meanwhile changes by orders of magnitude from its bare value. If g{A) flows to strong coupling, i.e. exceeds 
the cutoff Aq, at a scale As larger than the thermal scale At then there is a wide range of frequency scales 
over which the junction exhibits universal superconducting behavior. Similarly in the insulating side of (60), 
if 5(A) decreases by an order of magnitude at a scale A/ larger than the thermal scale then there is a wide 
frequency range over which the Junction exhibits insulating behavior. 

Let us now estimate the scales A/, Ag and At- To this end we first solve equation (60) to obtain g{l) = 

goexp (j^l — R)l^, where we have defined .R = (1 -|- F)R/Rq. When g is irrelevant (^ > 1) we get the 

typical (logarithmic) scale Ij = (R — 1)^^ and therefore Aj = Aoexp(— (i?— 1)~^). If on the other hand g 

is relevant {R < 1) it reaches the cutoff at the scale Ag = Ao(s(o/Ao)^^*-^~^^. Now to find the thermal scale 
we substitute g{l) into the flow equation (54) for the temperature. We then solve for T{1) with the initial 
condition of zero temperature T(0) = to obtain 

^ Ao 2^ - 1 V J Aq 2R-1 ' ^ ' 

The last equality is for large I and R > 1/2 which is the regime of interest in any case. Prom here we 
immediately obtain the thermal scale 



\2R-l) Ao 



Iff (62) 



The crossover between the superconducting regime and the thermal regime in the system parameter space 
R/Rq, F and 50/ Aq is defined by the surface At = As which is given by 



1-(1 + F)R/Rq 

go _ ( IqR/Rq \ 2(i+f-)R/RQ-i 



Ao \2{1 + F)R/Rq-1 



(63) 



At larger values of go the junction flows to strong coupling and displays superconducting behavior over a 
wide frequency range. 
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FIG. 2: Phase diagram of the noisy Josephson junction, as function of the bare resistance R/Rq and the bare 
Josephson coupling go/Ao, for different values of the noise: (a) F = 0.01, {b)F = 0.1, (c) F — 2.0. The blue lines 
(lower half of each subfigure) are obtained from the weak coupling analysis and correspond to equations (63) and 
(64). The red lines (upper half of each subfigure) are obtained through the duality transformation. The crosses (+) 
indicate the points used for the numerical solution of the RG equations shown in Fig. 3. 



We can find the crossover between the thermal regime and the insulating regime as the surface on which 
the thermal scale equals A/. This crossover is given by 

Aq IoR/Rq 

For smaller values of go the junction displays an insulating behavior controlled by the weak coupling critical 
point. For larger go, the effective temperature smears out the insulating behavior. 

If the bare Josephson coupling is large with respect to the cutoff scale Ag, one can apply the duality 
transformation described above and obtain an effective description of the junction in terms of a small phase- 
slip fugacity. In the dual picture, the roles of the insulating and superconducting regimes are exchanged. 
Therefore the dual version of eq. (63) describes the crossover between the insulating and thermal regime, 
and that of eq. (64) the crossover between the thermal and superconducting regimes. 

The phase diagram derived in this way is shown in Fig. 2 for both the weak coupling and strong coupling 
regimes for different values of the noise power. We see that the noise has two main effects on the phase 
diagram: first it shifts the transition from superconductor to insulator away from the universal value Rc — 
Rq. The shift is in opposite directions in the weak and strong coupling regimes. Second, increasing noise 
leads to a growing crossover region which exhibits effective thermal equilibrium behavior between the two 
states. In the next section we shall derive the universal transport characteristic expected in each of the three 
regimes. 

To confirm the phase diagram obtained from the approximate solution of the RG equations, we solve the 
equations (53-56) numerically for various values of the bare parameters. Figure (3) shows an example of 
the RG flow for the following set of parameters 50/Ao = 0.1, Ro/Rq = 0.7, Tq = and F = 0.01,0.1, 1.0. 
These points where chosen to demonstrate a noise-driven transition and are located respectively in the super- 
conducting, thermal and insulating regimes. We see that in the superconducting regime (a) the Josephson 
coupling grows rapidly and reaches the cutoff scale (A) before the temperature does. In the thermal regime 
(b) the temperature grows more rapidly than the Josephson coupling, preventing the latter from reaching 
the cutoff scale. Finally, in the insulating regime (c), the Josephson coupling decreases significantly, before 
the temperature becomes dominant. 



D. Current voltage characteristics of the junction 

The different dynamical regimes determined from the RG flow are characterized by distinct transport 
properties. Here we derive the I-V characteristics of the junction. We find a simple universal behavior in 
the different regimes that nevertheless betrays the non-equilibrium nature of the system. 
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FIG. 3: Example for RG flows of the Josephson coupling g{l) and of the effective temperature T{1) in the different 
dynamical regimes. We used the parameters (70/Ao ~ 0.3, R/Rq = 0.7, To — and different noise parameters (a) 
F — O.OOf , happens to be in the Superconducting regime. The Josephson coupling reaches the cutoff scale Ao and 
runs off to strong coupling before any appreciable temperature is generated, (b) F = 0.1 is in the thermal metal 
regime. The Josephson coupling may grow initially, but the effective temperature is first to reach the cutoff scale. 
At lower frequency scales the Josephson coupling is strongly suppressed by the effect of temperature, (c) F — 2.0 
is in the insulating regime. The Josephson coupling is irrelevant and it has a chance to decrease significantly from its 
bare value before the temperature reaches the cutoff scale. 



1. Insulating regime 

The insulating state occurs in the regime where the Josephson coupling is irrelevant. We therefore expect 
the transport properties to be well described by a perturbative analysis in the weak coupling limit. The 
standard expression for the current through a voltage-biased shunted Josephson junction, to second order in 
the Josephson coupling g is 

/ = /„ + = ^ + gHmUiuj = eV) (65) 

Here J„ and Is are respectively the normal current and the super-current through the junction. TZ{uj) is the 
response function of the operator e'*'-*^ in the steady state of the RC circuit unperturbed by the Josephson 
junction: it corresponds to the Fourier transform of (6). This perturbative calculation is sometimes referred 
to as the P{E) approach in the context of Josephson junctions"^^. For voltages significantly higher than the 
thermal scale A^, we can use the zero-temperature result (49) to obtain 

ls{V) = ^hV (66) 

Here Ii — Ii{R/ Rq, F) is a unitless constant, whose explicit expression is given in Eq. (A19) 

The above perturbative result can be improved by utilizing the RG approach. Instead of using perturbation 
theory on the bare model we apply it to the renormalized model at the scale A = ^ at which point the RG 
flow must be terminated. If we neglect the temperature, this gives exactly the expression (66), only with 
A V and renormalized values of parameters 

lAV) = hiWf-^ (67) 

Here Ii{l) — Ii{R{l)/ Rq, F{1)) and ly ~ log(A/V^). We see from this that the I — V curve is a power-law 
only at voltages higher than the thermal scale, at smaller voltages, the Josephson coupling is suppressed by 
thermal fluctuations and the contribution of cooper pair tunneling to the total current becomes negligible. 
In addition, the RG flow shows that the algebraic decay is actually controlled by the renormalized resistance 
R — R{1 — log(l^)) and noise F = F{1 — log(F)), rather than by their bare values. This correction is 
especially important if the bare Josephson coupling is large go ^ A and the insulating regime is reached only 
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FIG. 4: I-V characteristic of the Josephson junction in the insulating regime of the phase diagram Fig. 2, calculated 
from the RG flow using Eqs. (68). The bare parameters of the junction are the same as the one chosen for Fig. 3(c): 
(?/A = 0.3, R/ Rq — 0.7, and F = 2.0. If the junction was a perfect insulator all the current would flow through 
the resistor and the curve would be the dotted line V = RI. In reality the junction carries a non-vanishing current 
Is leading to a power law deviation from the linear curve. INSET: double logarithmic plot of the supercurrent 7^. 
The dotted line corresponds to an ideal power law with exponent R{It)/ Rq{1 + F{It))- The actual Is follows the 
power-law down to the effective temperature TeS = Ae~'^ 




FIG. 5: Sketch of the effective circuits used for calculating I-V characteristics, (a) is just the bare physical circuit. 
Panel (b) shows the renormalized circuit including the effective conductance 5G that is generated in the course of 
renormalization. Panel (c) shows the bare circuit corresponding to the dual theory and (d) is the renormalized dual 
circuit. Note that the dual current Iduai can be thought of as the current of flux quanta trasversing perpendicular to 
the junction and the resistor in panel (a). 

as a consequence of a long RG flow. If, on the other hand, the bare Josephson coupling is small go ^ A, 
both the corrections to R and to F are negligible, and the perturbative result coincides with the RG result. 

2. Superconducting regime 

In the superconducting regime the Josephson coupling runs off to strong coupling well before reaching the 
generated thermal scale. Therefore a perturbative analysis of the superconducting regime can be made from 
the dual standpoint in which the phase slip fugacity is the small parameter. 

We start from a direct perturbative analysis applied to the bare dual circuit. To understand the connec- 
tivity of the dual circuit, note that a voltage bias of the original circuit corresponds to a constant current 
^duai = V/e of flux (dual particles) which flows perpendicular to the original circuit, first across the junction 
and then across the resistor. The current through the original circuit corresponds, on the other hand, to the 
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FIG. 6: I-V characteristic of a noisy Josephson junction in the superconducting regime for (?duai/Ao = 0.3, Rq/R = 
0.7, F = 2.0. In the regime Te// <^ 7/2e <^ Ag,duai we see the algebraic behavior, predicted by perturbation theory 
in the strong coupling limit (69). In the regime 7/2e <C Teff , highlighted in the inset, the junction follows Ohm's law 
with an effective resistance 5R — SGduai , perturbatively smaller than the bare resistance R (dashed line) 



sum of the dual voltages on the junction and on the resistor. Hence, in the dual picture we have a resistor 
-Rduai = Rq/R connected in series with a Josephson junction gduai <C Aq. This is sketched in Fig. 5. (See 
also appendix B for a detailed derivation of this result). 

We can now apply the weak coupling result (66) to obtain /dual versus T4.duai on the dual Josephson element 
alone. Then using /dual = V and l^^.duai = / and substituting the dual junction and noise parameters we 
can translate the result to physical voltage versus cooper pair current on the original Josephson junction. 



^ /? ' y A V A , 

Before applying the RG approach to compute I-V relations we should understand how the connectivity 
of the dual circuit changes in the course of renormalization. Under the renormalization we generate an 
effective dissipation, which can be regarded as added conductance in the circuit. The (dual) current that 
went through the bare Josephson junction, now splits between the junction and the conductance added in 
the renormalization process. This implies that the added effective dissipation is connected in parallel to 
the dual Josephson junction and both are connected in series to the resistor Rduai{0) of the bare circuit as 
depicted in Fig. 5 (b). 

Now to compute the I — V relation we take the voltage on the dual circuit as the cutoff to the RG flow 
A* = Vduai ■ The total (dual) current is then given by: 

Idual — SGdual{l*)Vdual + Il{h) ^'^'t'f^ * (70) 

Vdual 

Where 5G{1) is the renormalized conductance in parallel to the JJ at the scale I. The last relation can be 
immediately translated to one between the physical voltage and current 

V^SGdualiQl + hiU f'^ff (71) 

Note that SGduai is proportional to Iidduai/ hence both contributions to the voltage on the Josephson 
junction are of order gd^ai might have been expected. The first term in (71) constitutes a small linear 
resistance, that accompanies the emergence of effective temperature and dominates the transport at very low 
currents (//2e < T^jf)- In the other limit, when the current is much higher than the effective temperature, 
the second term in (71) dominates and leads to essentially the same result as bare perturbation theory (69). 



E. Impurity in a Luttinger liquid 



As we have seen in Sec. II B, the non-equilibrium action of a local impurity in a one dimensional system 
driven by 1/f noise can be mapped exactly to the action considered here (2). This mapping allows us to use 
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the results from the previous section to determine the non-equihbrium phase diagram and IV curves of the 
one dimensional case. The same mapping allows us to compute the current-voltage characteristics as well. In 
the Id model, it is possible to induce a constant flow of particles by the transformation (j){x, t) — ?■ 0(x, t)+TTlt. 
(Recall^° that the displacement field </'(a;o, t) jumps by n each time that a particles crosses the point x = xq). 
The flow of particles (passing through both the one dimensional system and the impurity) is therefore 
equivalent to a voltage over the shunted Josephson junction (acting on both the junction and the transistor). 
To introduce a constant voltage in the Id system, we need to add the following term to the Hamiltonian 

V j dx p{xo) = ^2</)(a;o) (72) 

Here p{x) « {—l/(j))dx(t>{x) is the density of particles at position x. Eq. 72 shows that the voltage of the 
Id system (given by the sum of the voltage over the Id lattice and the voltage over the single impurity) 
corresponds to the current through the noisy Josephson circuit (given by the sum of the current through the 
resistor and the super-current through the junction). Using this mapping, we find that current- voltage (IV) 
characteristic of the noisy Id system with a local impurity is precisely identical to the voltage-current (VI) 
curve of a noisy shunted Josephson junction that we computed in Sec. IV D. 

V. ONE DIMENSIONAL SYSTEM 

We turn to discuss the model of interacting particles in one dimension in a commensurate periodic potential 

subject to 1/f noise. The model was introduced in Ref. 8 and reviewed in Sec. II B. We have pointed out 
that in absence of the periodic potential the system exhibits algebraic crystalline correlations and response 
functions that may be very different from the equilibrium ones. The power law decay of the correlations 
depends on the dimensionless ratio between the noise strength and the coupling to the bath, F = ■7r~'^F/ {urj). 
Recall that we are considering the limit where both r] and F tend to zero but the ratio between them F may 
be arbitrary. 

Having found a scale invariant state in the quadratic model, we should ask how the non-linear coupling of 
(8) behaves under the same scaling transformation. In Ref. 8 we read off the anomalous scaling dimension 

of the cosine term from the power-law decay exponent of the crystalline correlations (11), which gives 
dg/dl = g {2 — K {l+ F^Y This result suggests a possible noise-tuned phase transition in which the periodic 
potential g changes from being irrelevant to being relevant. In equilibrium [F = 0) such a quantum pinning 
transition indeed occurs, and as can be seen from the scaling of g, it occurs at the universal value of K = 2. 

The purpose of this article is to go beyond the simple scaling argument and account for the back action 
of the periodic potential on the steady state correlations. This will be done using the general RG framework 
laid out in section III. As mentioned already, there is a formal similarity between this problem and that of 
the steady state following a sudden change of the Luttinger parameter in the Sine-Gordon model. The RG 
equations for the latter problem were derived in Refs. 9,10. We shall briefly rederive these equations here 
for completeness and analyze the dynamical regimes they give rise to in the noisy Luttinger liquid. 

A. RG equations 

The general framework that has been laid out in Sec. IIIB can be directly applied to the one-dimensional 
system. First we derive from the general formula (26) an explicit scaling equation for the coupling g valid 
to first order in the coupling. Then we use the general formula (40) for the scaling to second order in g. We 
carry out an expansion of this result in small frequency u and wavevector q in order to recast the functional 
equation (40) as scaling equations for the set of parameters in the action (8). Readers uninterested in the 
technical details of the derivation can skip directly to equations (82)- (87) that summarize the result. 

One technical note we should make before proceeding with the, is that for the derivation in section III B we 
assumed a cosine coupling of the form cos{(j){x , t)) , while we are now interested in the coupling cos(20(x, i)) 
appearing in the action (8). We should therefore make the substitution (f) ^ 2(j) in the formulas appearing 
in Sec. IIIB when considering this problem. 

To find the first order flow of g explicitly, we need to compute the variation of the fluctuation {(p) within 
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the quadratic theory changing cutoff scale. Neglecting the effect of the dissipation r] we get 

OAv^Oj |_1_(^2 _y2^2) ^^^^1^ 

= if /(t. dq) dLO u\q\ e-"l«l/^ ( cotgh + d{uj' - {uqf) (74) 

(75) 



if 



T 

l + -r+F 

A 



Here we used an exponential cutoff function Tx{q,uj) = cxp{—u\q\ / A) and F is defined in (10). Substituting 
this result in (26) leads to the scaling equation given in (82) below. 

We now turn to derive the flow equations to second order in g through an expansion of the functional 
equation (40) in small q and lu. A contribution with both q = and a; = is found only in the Keldysh 
component dGj} {q, ui) and is given by 

^2 f f „2 



'A2 



i2-2Dg) Jd{x/u) J dtK:{x,t) = -2i^lT{F,K) (77) 

Here Dg is the scaling dimension of the commensurate potential found from the first order flow for vanishing 
temperature and dissipation, Dg = 2 — K{1 + F). JC{x,t) is the correlation function of the operator e*^"^ 
taken in the quadratic action with T = and 77 = 0. We used the fact that under these conditions the 
correlation is a scaling function of Aut and Ax, K = IC{Ax/u, At). As discussed already below Eq. (48), this 
contribution to the Keldysh component can be associated with an effective temperature in the low frequency 
limit through G^^(0, 0) -H- rjT. The identification is made by comparing to the low frequency expansion of 
the Keldysh component at equilibrium r]UJ coth(a;/2r) = rjT + 0{uj'^). 

By inversion and time reversal symmetry dGj}{x,t) does not contain terms proportional to q or oj, and 
its next non- vanishing terms are proportional to q^ and ui^. These terms are less relevant than the bare 
quadratic action, and therefore they do not affect the correlation and response functions at long times. 

In contrast to the Keldysh component, the renormalization of the off-diagonal component of the action 
dGji{q,u)) does not have a zero frequency component. The lowest order term in the gradient expansion of 
this term in (40) is 

- ^^Y^ {^Jl^'t^")-^) jdttn{x,t) (78) 

{l-2Dg) j d{x/u) J dtn{x,t) = iijj^In{F,K) (79) 



'A2 



This contribution amounts to renormalization of the dissipation as drj/dl = ig^/A'^I^{F,K). 

The next order in the gradient expansion of dGj^^ gives terms proportional to q'^ and w^. These terms 
renormalize the bare term {uK)u)'^ — {u/K)q^). The independent renormalization of K and u are then given 

by 

dK~^ f f - 

__ = -2dlj^^Dg j d{x/u) J dt {{x/uf + t^mx, t) ^ dl ^Jk{F, K) (80) 



^ = -2dl^^Dg j d{x/u) j dt {{x/uf - t^)n{x,t) ^ dl -l^I^{F,K) (81) 
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The above results are summarized in the following RG scaling equations: 



dl 

djyF) 

dl 
djyT) 

dl 
dK^ 
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drj 



dl 
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Ao ^ TT^v) 



T F \ 



(82) 



(83) 



du 
Hi 



■''■0 

2tjT+^It 



(84) 



(85) 



(86) 



(87) 



The cocfRcients It, Iri, Ik, Iu arc defined in (77-8f). Analytic formulas for these coefficients as ftmctions 
of K and F are derived in Appendix A 2. It should be noted that precisely the same scaling equations were 
obtained in Ref. 10 in the context of a quench of a one dimensional system in the presence of a commensurate 
potential. (See also Sec. II B). 

Let us briefly point to the key properties of the above scaling equations. The first order flow of the coupling 
g described by Eq. (82) is very similar to the equilibrium situation but with a negative contribution of the 
noise through F to the anomalous scaling dimension. It is important to note that at this order in the RG 
the noise, although weakening the correlation functions, is not acting as a temperature, since it preserves the 
power-law behavior of the correlations. One might thus expect some change in the dynamics of the particles 
from delocalized to localized when the noise is tuned to change g from being irrelevant to relevant. At this 
level, the physics is identical to that described by the simple scaling argument in Ref. 8. 

As in the case of the Josephson junction new physics is introduced by effects that are second order in the 
coupling g. Eq. 83 for example leads to the emergence of finite dissipation rye// ~ Iri9^ I in the steady state 
although the bare dissipation (r?(0)) was infinitesimal. The analytic expression given in Eq. A29 shows that 
the coefficient vanishes identically at equilibrium F = 0, but also that in presence of the noise the effective 
dissipation diverges for K = K{1 -\- F) < 2, that is exactly where g becomes relevant. We can interpret this 
result as emergence of sub-ohmic dissipation that is not captured by the perturbative RG scheme. In what 
follows we restrict ourselves to the regime K > 2, where rjes can be made as small as we want by taking a 
sufficiently weak lattice. The renormalization of the dissipation generates also an effective renormalization 
of the noise (84), introduced to conserve constant the Keldysh part of the action. 

The emergence of dissipation is accompanied by the generation of an effective temperature described by 
Eq. (85). We can estimate this temperature from the effective fluctuation-dissipation relation between the 
Keldysh and retarded parts of the action 



This is a rather suprising result as it implies that the temperature generated in the steady state may be 
large even for arbitrarily small non-linear coupling g and noise F . Furthermore, if we track the generated 
temperature as a function of the flow parameter we see that the effective temperature reaches its large 
value immediately at the beginning of the flow. That is, after one inflnitesimal RG iteration we have 
T(Z = 0+) = Toff. From this point of view it looks like the non-equilibrium critical state that exists in the 
purely quadratic model g — is an extremely singular limit. We shall discuss below the implications of 
this result and argue that the critical state of the quadratic model does nevertheless in practice govern the 
physics of the system with a weak non-linearity. 

Note that in the case of the Josephson Junction discussed earlier the effective temperature was paramet- 
rically small in both the noise strength and the non-linear coupling. The difference between the two cases 
stems from the different nature of the dissipation. While the effective dissipation in the one-dimensional 
system is small, of order g^, in the Josephson junction model even the bare dissipation is of order 1. 

The scaling equations (86) and (87) for the dimensionless coiiplings K and u are controlled by the integrals 
Ik and /„, defined in (A30) and (A31). Though we do not have analytical expressions for these coefficients, 
numerical evaluation of the integrals shows that both are finite for any K > 2. Moreover the factors Ik 



Teff « IrgViVeff^^) = Ao/t//^ ~ Ao(^ - 2) 



(88) 
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and lu do not vanish in the equihbrium case F = 0. This is expected given that K should also have a non 
trivial flow near the Kosterlitz-Thouless transition of the Sine-Gordon model at equilibrium (see for example 
Ref. 19). 

B. Non-equilibrium phase diagram 

Analysis of the RG equations (82-87) suggests that there are at least two distinct dynamical regimes in 
presence of 1// noise {F > 0). First, for K = K{l + F) < 2 the dissipation 77 diverges through the divergence 
of the coefficient This probably indicates generation of local sub-ohmic dissipation which is beyond the 
present RG scheme. For ^ > 2, an effective temperature T^s ^ {K — 2)Ao is generated immediately at 
the outset of the flow. We will now consider each regime separately and then move to the crossover region 
between them. 



1. Critical - thermal regime 

To understand the nature of the thermal regime, recall that the; source of thermalization in the system is 
the mode-coupling term g, which scatters the non-equilibrium population of high-frequency modes towards 
lower frequencies. The rate by which this scattering occurs is the generated dissipation rate ryeff « I^g^ j 
This rate sets the time scale for the system to achieve the steady state if the non-linearity were turned on at 
some point in time. But, once the system reaches steady state, its effective temperature Tgjj depends only 
on the energy density in the initial critical steady-state and not on g. Thus the singularity of the state at 
(7 = is resolved by the time-span over which the experiment takes place. In order to see the true steady 
state with temperature T ~ {K — 2)Ao we would need to conduct an experiment over a time that diverges 
as T ~ l/g^ as g — >■ 0. As a consequence, if we restrict to time scales shorter than l/r/cfr, the system is 
governed by the quadratic part of the action (8) with negligible T. This is a quasi steady state, with scale 
invariant correlations and response given by (11) and (12). 

Z. Local dissipation 

As the effective Luttinger parameter decreases and approaches K = 2, the renormalization of the dissipa- 
tion diverges and the dissipation constant inevitably becomes the first to reach the cutoff. The scale at which 
this happens is the generated effective dissipation. Clearly at that scale the dissipation is the largest scale 
in the system and we can neglect the original kinetic term K'^^co'^ compared to it. Therefore the effective 
action we should work with, below the dissipations scale A^ = rj^s is: 

This describes the thermal regime of a critical theory with dynamical critical exponent z = 2. It is interesting 
to note that in this theory the noise becomes a relevant perturbation with scaling dimension Dp = 1, exactly 
like the temperature. 

3. Crossover regime 

We have identified above two dynamical regimes. For K < 2 the renormalization of the dissipation 
diverges, which leads to a thermal state with short range spatial correlations already at initial stages of the 
flow, that is even at frequency scales. On the other hand, ior K ^ 2 there is a wide range of frequency scales 
in which the dissipation is still small. If the finite time of the experiment does not allow for equilibration to 
the high effective temperature of the steady state then the physics will be governed by the non equilibrium 
critical state of the Luttinger liquid. 

Of course at sufficiently small frequencies, the generated dissipation being a relevant perturbation will 
eventually dominate. So, the transition between the two regimes is not perfectly sharp but rather a smooth 
crossover. Nevertheless we can chart the line of the crossover in the space of K versus the non-linearity g by 
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FIG. 7: Crossover line (93) as function of the Luttinger parameter K and the commensurate lattice strength g/A^, for 
different values of the noise strength ■K~'^F/ri: (a) 0.01; (b) 0.1; (c) 1.0. In the "local dissipative" regime the effective 
dissipation is diverging and the effective temperature is small, while in the "critical thermal" regime the effective 
temperature is large and the dissipation small. A critical behavior is still observable in the short-times dynamics 

invoking similar arguments we used above for the Josephson junction. Namely, we should find the regime in 
K, g space, where a critical flow of g down by an order of magnitude can be observed before the dissipation 
reaches the cutoff scale, where g changes by an order of magnitude. 

The typical scale at which g is reduced by an order of magnitude in absence of dissipation is found by 
solving the scaling equation for g, assuming no renormalization of the dissipation zero temperature and. The 
solution 



5(0 = 5oe 



(90) 



, gives the frequency scale Ag = Ae*^^ . Now to determine the typical scale where dissipation enters 
significantly we solve the scaling equation for r/ (85) while substituting the above solution for g. This gives 



77(0 = Iv 



A3 



„2(2-/f)/ 



(91) 



" 2A3 - 2 



(92) 



Here we assumed 1^ to be constant, neglecting the renormalization of the Luttinger parameter and of the 
sound velocity. Recalling that 1.^ cx {K — 2)~^, we find that as we move away from the transition point 
K = 2, the effective dissipation A^ decays like {K — 2)~^. 

Comparing A^ and Ag, we obtain an analytic expression for the crossover line between the dissipative 
regime and the critical regime: 



g2 2if-4 



A4 



(93) 



For smaller values of g, the dissipation scale is smaller than the scale for g in the critical state and it is 
possible to observe the critical regime. At i^T —>■ 2, the transition line goes to zero: any infinitesimal amount 
of non-linearity is enough to generate an extremely strong dissipation. At large K the effective dissipation 
is small, but, as will see below, the effective temperature becomes large. 



VI. SUMMARY AND DISCUSSION 



In this paper we investigated the physics of low dimensional systems driven out of equilibrium by external 
noise sources. We considered two specific models which in equilibrium undergo interesting phase transitions: 
(1) a resistively shunted Josephson junction and (2) a one-dimensional quantum liquid in a commensurate 
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lattice potential. The basic question underlying our study is whether non equilibrium systems can exhibit 
universal physics. We have recently shown* that such systems, in the absence of non-linear coupling between 
the modes, establish a critical steady state when driven by 1// noise. Here we used this scale invariant 
steady state as a starting point for a controlled renormalization group analysis to treat the effect of the non 
linear couplings and thereby identify universal aspects of the non equilibrium behavior. The perturbative 
RG scheme is formulated within the Keldysh path integral framework. 

The scaling equation at the first order in the coupling constant indicates a transition between different 
regimes where the scaling dimension of a relevant operator controlling the physics of the system (Josephson 
coupling for the noisy Josephson junction, periodic potential of the lattice for the ID quantum liquid) 
changes sign. This occurs as a direct extension of the quantum phase transitions that would exist in the 
corresponding system in equilibrium. Here however the scaling dimension can be tuned by varying the 
external noise strength and not only intrinsic parameters (i.e. shunt resistance or Luttingcr parameter). At 
this level the RG analysis is consistent with the scaling argument we gave in Ref. 8. It is important to note 
that at this order, although the noise is weakening the various correlation functions of the system, it is acting 
in a very different way than a temperature, and in particular preserves the power-law decay of the various 
correlations in the system. 

The back action of the nonlinear coupling on the steady state correlations is only given by the RG 
equations at the second order in the coupling. We find that generically this leads to the emergence of an 
effective temperature and effective dissipation in the steady state. In the case of the Josephson junction 
both effects arc pcrturbatively small, allowing the observation of a universal crossover between an insulating 
and superconducting regimes, controlled by the above non-equilibrium quantum critical state. Using the RG 
approach, we obtain analytic expressions for the shape of the crossover, summarized graphically in Fig. 2. 
The RG approach allows us to compute the physical properties of each regime, such as the IV characteristic 
of the junction shown in figures 4 and 6. In general, the IV curves show an algebraic behavior, determined 
by the underlying critical state, followed by a thermal Ohmic behavior, at voltages lower than the effective 
temperature. 

In the one dimensional case, we considered two types of perturbations around the critical state. The first 
type, a local impurity, can be mapped into the noisy shunted Josephson junction described above: the one 
dimensional system plays the role of the resistor and the impurity plays the role of the Josephson coupling. 
The second typo of perturbation, a commensurate lattice, leads to new physics. Despite the different initial 
set-up, the RG flow is equivalent to the one dimensional quench considered in Ref. 9. The resulting non- 
equilibrium phase diagram shows two distinct regime: a thermal regime and a dissipative regime. In the 
thermal regime, the lattice is irrelevant at first order. However at second order in g the coupling of modes 
at all scales induced by the lattice acts to generate fluctuations and effective dissipation, both proportional 
to ff^ . This implies a fluctuation dissipation relation with an effective temperature of the order of the cutoff, 
from which the coupling g drops out. This seemingly paradoxical result is coming from an inversion of two 
limits, namely the fact that the coupling to the noise can be small and the fact that one observes the steady 
state of the system at infinitely long time, allowing a large amount of energy to be pumped into the system 
under the form of a temperature. However, as the effective dissipation is small, the thermalization rate of the 
system is vary small, thus allowing to observe transient non-equilibrium phenomena for relatively long times. 
Clearly, the complete study of the dynamics of the system, in particular at transient times goes beyond the 
scope of this paper and of the perturbative RG analysis, and deserves further study. 
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Appendix A: Explicit calculation of the RG flow coefficients 

In this appendix, we show how to define a proper cutoff and use it to compute the dimensionless pre-factors 
appearing in the scaling equations (54-56) and (83-87). Specifically, we need to obtain explicit expressions 
for the correlation JC{x, t) and response 7?.(a;, t) of the exponential operator e*'^ and then use these expressions 
in order to compute the integrals in (48-51) and (79) and (81) 
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FIG. 8: Real and imaginary part of the cutoff function (Al). The real part of this function equals to e '"1/^. In 
addition, the cutoff function has a non trivial imaginary part which allows to preserve causality. 



1. Noisy Josephson junction 

The simplest choice for the cutoff would be a step function in frequency, r(a;/A) = 1 — 9(w/A). However, 
as we noted earlier, this choice is not valid because it violates the causality of the retarded Green's function 
Gr.a{^) = Gij^Ao ('^)rA(w). Instead we choose the following cutoff function 



-\u^\/A 



{v - bj - iQ+){v + uj + 

2i [ ^'^\ ( '-^ 

Coshlntegral I — ) Sinh I -r ) + Sinhlntegral I — ) Cosh I — 



(Al) 



Eq. (Al) explicitly preserves causality and, as shown in Fig. 8, satisfies the defining conditions of a cutoff 
rA(0) = 1 and r((X)) = 0. 

This specific choice allows us to obtain simple expressions for the Green's functions, GK,A{t) = {(l>{t)ip{0)) 
and Gfl,^A(t) — {(j){Q)(l){t)) . For convenience, we first compute the derivatives of these functions with respect 
to A. Assuming zero temperature (T = 0) we obtain: 
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Here Q{t) is the Heaviside step function. We can now integrate with respect to A and obtain: 



(A2) 
(A3) 
(A4) 
(A5) 



GK,Ait) = GK,A(0) + 2i-^(l + F)log(l + A2t2) 

R ( t 

GR^Ait) = 2e(0-^arctan(^- 



(A6) 
(A7) 



The correlation and response functions of the exponential operator e*^ are in turn simple functions of the 
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Green's functions Gk,a and Gr^l'- 

j^(^) = ^gie(t)+i9(t) ^g-i9(o)-j0(o) ^gie(o)-je(o)^^ ^^g-j 

= g-i((o(o+e(t)-e(o)-e"(o))") _^g-i((e(t)+e(t)-e(o)-e(o))'> ^^9^ 

= 2cos(Gj{,A(i'))e'^^^'^^*^"°^'^^°^^ (AlO) 

= ^gie(t)+i0(t) ^g-ie(o)-ie(o) _ gie(o)-ie(o)^^ j-^-^^^ 

= 2sin(Gfl,A(t'))e'^^''"'^*^"^'''''^°'' (A12) 

Here we used Gk = i{O0), Gr = i{ee), and {66) = 0. 

We now use these expressions together with the first order result (41) to derive the Iq and I\, defined in 
(48) and (51) 

/. = (2 A„ + _ 1) £ (2^„., (A,)) ( (A13) 



2 — (l + F)-ljy ^^dy cos(^2— y) (cosy)'«Q^'+"^ ' (A14) 
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(A15) 



7. = (2^(1 + ^) - 2) £ * sin (2 A,,tg (At)) (^) (A16) 
= 2 (^2-^(1 + F)- 2) ^"^'dy sin(2-^y) (cosy)'^^'+^^-' tan(y) (A17) 
4^ 1"'^ dy cosf2^y) (cos y)'^^'+-^'"^ (A18) 

(A19) 
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Here B{x,y) = T{x + y)/T{x)T{y) is the Bigamma function. This function diverges if a; = and y 7^ 0, 
implying that both /q and Ii vanish at equilibrium {F = 0). In addition, we note that eq. (A19) is valid 
only for 2;^(1 + F) > 1, while for 2;^(1 + F) < 1 the integral diverges. 

2. One dimensional system 

We now move to the one dimensional system. Here the natural choice of the cutoff is r(g,w) = e""'*'/^. 
The correlation and response function of ^ are then given by: 

GK{x,t) - GK{0,0)=i^mx,t)-cl){0,0)f} 



= i^Jiu dq) da; (1 - e^"*-"'"^*) e""!'!/^ (^1 + -^ij) 5(^2 _ ^2^2^ 
= (1 + F) [log(|l+A(x/u-f)|) + log(|l+A(x/u + i)|)] (A20) 

Gnix^t) = iU(x,tU(0,0)) = lim / (u dq) duj e'«^-'"* e-"l«l/^^5 (A21) 

r)->o J u-^ — u-^q-^ — irjUJ 

= y [artg(A(x/w + 1)) + artg(A(a;/u - t))] (A22) 

Here we have neglected the effects of both the dissipation and the temperature (r/ = T = 0+) and F is 
defined in (10) 
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Substituting these expressions in (77-81), we find: 

/CO POO / 1 \ / 1 \ ^^'^ 

dx dt COS {K aitan{x + 1) - KaTtan{x-t)) [- ( :j r^] (A23) 
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In the absence of the noise [F/r] = => K ^ K), both It and go to zero. As expected at equihbrium, 
the non-hnear couphng does not lead to a renormahzation of the temperature or of the dissipation. Also 
note that, for any finite noise strength, 7^ diverges at ^ = 2, while 7^ remains finite at this point. This 
observation will have dramatical consequences on the properties of the expected non-equilibrium phase 
transition^, located precisely at TsT = 2. 

For the integrals Ik and 7„ we were not able to find an analytic solution. However, we have numerically 
computed them and found that, in agreement with the known equilibrium results^^, both integrals do not 
vanish even in the absence of the noise. In addition, we found that both Ik and 7„ are finite for any K > 2. 
For details about their behavior around the K = 2 point, the reader is referred to Ref.lO. 



Appendix B: Duality transfortmation 

In this appendix we review how the duality transformation''^'''^'''''' works within the real time Keldysh 
framework. Our starting point is the non-equilibrium Keldysh action (2), which we split into a unitary part 
and a dissipative or non-local part S = + Sdis- The unitary part of the action includes the capacitor C 
and the Josephson coupling while the non local terms include the resistor and the external noise. 

In the strong coupling limit {g ^ A) the unitary action is large justifying the use of a Saddle point 
approximation for the phase dynamics of the field 9. This leads to two independent equations of motion for 
the fields ^± (t) on the forward and backward branches of the Keldysh contour 

^ = {RC)d^Mt) - 5sin(0±(i)) = (Bl) 

Equation (Bl) has instanton solutions (j>{t) = 2Tresf{t — to), where = ±1 and f{t) is a step function 
broadened by a short time-scale r = \fcjg. The instanton solution describes a jump of the phase (phase 
slip) by ±27r between minima of the cosine potential. A more general saddle point solution can include any 
number n+ of instantons in the forward and n_ instantons in the backward branches of the Keldysh contour. 
Furthermore, in the strong coupling limit phase slips are rare, in the sense that the average time between 
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them is much larger than the instanton time t. In this case we can take a general saddle point solution to 
be a sum of independent instantons and at the same time regard each instanton as a step function 



n± 



</'±(i) = I^e±,ie>(*-i±,i) (B2) 

Here t±^i is the time of the instanton. Plugging in these solutions into the Keldysh action, the path integral 
is converted to a sum over all possible instanton configurations on the two branches 

/ V<f>+Vct>- e'S[*+'^-l = E E A (^«"«5i„s + iSdis[c+, c_]) . (B3) 

■' a=± c. 



Here c± = {n±, e±,i, t±,j} denotes an instanton configuaration on the forward/backward branch and Si^s is 
the action of a single instanton. 

The non-local part of the action S'dis[c+, c_] can be written explicitly in terms of the instanton coordinates 
ti^a and Ci^a = i by plugging the instanton solutions into the the quadratic form of the original action action 



(B4) 

This part of the action contributes a pairwise interaction between instantons at different times. 
We now introduce the dual field 9 ait) through the identity 

eiSai.[..,t.] ^ fmexpUY.apfduj u;2r(-w)GdisHa/3^''(w) + ^2^ae„,i0«(fi) ] (B5) 

Substituting this back into (B3) we can perform the summation over as follows 

E + ^ E ^Mk)^ = E (E e^^-'+^^^^'W^ (B6) 

= exp e'^-' cos(^„(i))^ = exp (^i ^duai J dt (cos(e„(t)))^ (B7) 

Here we used the time scale RC to transform the sum over times into an integral. The phase slip fugacity 
5duai = (l/i?C)e*'^=°' be estimated using the WKB approximation^^ to be (/dual ~ {l/RC)e^^^^'^\ 
Finally, putting (B7) together with (B5), we obtain the dual action 

S = J2[doJ ea{-u;){G^l,r^9p{u;) + <?duai / dt cos{e+) - cos(0_) (B8) 

We have defined G'^^^i{uj)ai3 = ct/3 Gdis(w)a,/3- 

To compare the dual model (B8) with the original one (2), it is useful to express G^Jj^i in terms of its 
classical and quantum components: 



<^dual <^dual / V '^Im ^I: 



Using (2) we obtain: 



/ i^co \ 

Cduai = (^a;cotgh(w/2T) + i^w^-aj J (BIO) 

By inspection we see that the dual action is identical to the original one, up to the duality transformation 
R/Rq Rq/R, F [R/ RqYF , T ^ T, and g S'duai- In this sense the shunted Josephson junction is 
self dual even in the presence of the external noise. 
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In the above derivation, wc implicitly assumed that the junction is disconnected from external sources 
apart from the noise source. However to derive the current voltage characteristics we need to embed the 
junction in a circuit with applied current or voltage bias and to generalize the duality transformation to this 
case. The question is how the circuit structure transforms under the duality transformation. 

Consider the circuit shown in Figure 5(a), in which the noise-driven shunted Josephson junction is biased 
by a DC current. The current can be enforced by adding the term Si — I J dt (t>{t) to the Keldysh action (2), 
which transforms the periodic cosine potential into a wash-board potential. If the current is small compared 
to g, we can still substitute the soliton solutions obtained above into the action. Doing this for the new term 
Si we obtain Si = 2TrIJ2a=± Si oi^iaUa- Introducing the dual field as before and summing over n, we get 

S = Y,jd^ w'^a(-w)(Gd4i)a,,3e/3(w) +9 j dtcos{9+it) + RqR) - cos{e-{t) + RqR) (Bll) 

(here we work in the units for which h = 1 and Rq = h/{2e)'^ = 2it). Equation (Bll) shows that a current 
I acts in the dual circuit as a voltage bias Vduai = RqIi applied only to the Josephson junction (and not 
to the resistor). In this sense, the dual model corresponds to a circuit in which the Josephson junction is 
connected in series to the resistor (See Figure 5 (c)). 

We will now show that the duality transformation maps the voltage across the original circuit into the 
current of the dual circuit. The voltage can be obtained formally by taking the derivative of the Keldysh path 
integral with respect to separate current sources 7+ and /_ acting in the forward and backward directions 
and setting /+=/_= 7 in the end: 

= Rq 5duai5t j d9 [sm{e+ + RqR) + sin(e_ + RqR)] t e'^'^-' (B12) 

As the junction is in a steady state, {V{t)) is independent on t and we can take for example t = 0: 

{V) = Rq 5duai j d0 [sm{0+) + sin(e_)] e*^^-' = Rq{ILJ (B13) 

We see from the last expression that {V{t)) = 7?q7|^jjj by definition of the super-current trough the dual 
junction. Note that since the Junction is in steady state the expectation value cannot be time dependent 
and we might as well set t = in Eq. (B13). 
Equations (Bll) and (B13) can be used to evaluate the IV characterstic of the original circuit, in its strong 

coupling regime. (The calculation follows the same steps used to compute the current across the junction in 
the weak coupling regime.) In the strong coupling regime g ^ A ffduai A and we can apply the first 
order perturbation theory in gduai to find 

(V) = i 5duai f dt'{[cos{e+{t') + RqR') - cos{e_{t') + RqR')] [sin(0+(O)) + sm{0_{Q))])^ (B14) 

J — CO 

Here (...)q is the average with respect to the quadratic part of the Keldysh action (B8). Finally, we exploit 
causality to extend the integral from t' = to t' = oo and re- write (B14) as 



{y) = .9duailni 



(B15) 



.J — OO ^_ 

^ gi^,M[n^ = RQl)] (B16) 
Here TZ{oj) is the Fourier transform of the response function of the operator e*". 
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